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ABSTRACT 

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of uniden¬ 
tified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing 
system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remain¬ 
ing unidentified gamma-ray source without a known association in the second Fermi-LAI source catalog 
(2LGL) and was among the top ten most significant unassociated sources in the recent third catalog (3LGL). 
PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of 8.9 Hz, a characteristic age 
of 49kyr, and spin-down power 1.0 x 10 36 erg s _1 . In 2009 August it suffered one of the largest glitches de¬ 
tected from a gamma-ray pulsar (A/// « 4.5 x 10 -6 ). Remaining undetected in dedicated radio follow-up 
observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location 
of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the 
interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects 
which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by 
which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for 
estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this. 

Subject headings: gamma rays: stars — pulsars: individual (PSR J1906+0722) 
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1. INTRODUCTION 


The large collecting area and continuous observati on m ode 
of the Fermi Large Area Telescope (LAT; [Atwood et al.|2009 ) 
make it an ideal instrument for the detection and analysis 
of periodic gamma-ray emission from pulsars. Through the 
careful analysis of the arrival times of photons covering the 6 
years since its launch, the LAT has discovered pulsed gamma 


ray emission from more than 160 pulsars 64 (Abdo et al. 2013 
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Carave o|2014[ ). 

While the majority of these pulsars wer e first found in ra¬ 
dio observations (e.g., Abdo et al.|2009b|c] ), the ephemerides 
from which could be used to test for gamma-ray pulsations, a 
substantial fraction of the gamma-ray pulsar population was 
discovered through blind searches of Fermi -LAI data (e.g., 
[Abdo et al.[2009at|Saz Parkinson et al.|2010|) 

In a recent work (Pletsch & Clark 2014| ) we presented 
newly advanced methods designed to increase the sensitiv¬ 
ity of blind searches without increasing the computational 
cost. These improvements have since been incorporated into 
a new blind survey of unidentified, pulsar-like Fermi-LAI 
sources being conducted on the distributed volunteer com¬ 
puting system, Einstein @ Home [^] Previous surveys have 
been extremely successful in det ecting new gamma-ray pul¬ 
sars ([Pletsch et al.] |2012a|b[ |2013| ), and the newly improved 
search methods, in combination with the latest Fermi-LAI 
data, offer a significant increase in sensitivity. 

As part of this survey, we carried out a blind search 
for pulsed emission from a point s ource in the third 
Fermi-LAI source catalog (3FGL, Acero et al. 2015| ), 
3FGL J1906.6+0720. This source, previously known as 
2FGL J1906.5+0720 ( [Nolan et al.||2012|) , is highly signif¬ 
icant and stands out as the most significant unassociated 
2FGL source. Moreover, it was included in the “bright” 
pulsar-like source list described by Romani ( 2012| . An in¬ 
vestigation of the spectral properties of 2FGL sources found 
that, after the source associated with the Galactic Center, 
2FGL J1906.5+0720 was the unidentified source most likely 
to contain a pulsar ( Lee et al. 2012 ). As such, over recent 
years, this source has been searched for pulsations, both in 


gamma rays (e.g., Pletsch et al. 2012a ||Xing & Wang|2014 ) 
and in radio observations (e.g., parr et al.| |2013j ). However, 
despite these attempts, pulsed emission from this source re¬ 
mained undetected until now. 

Here, we present the discovery and follow-up study of 
PSR J1906+0722, a young isolated gamma-ray pulsar de¬ 
tected by the Einstein@Home survey. 


2. DISCOVERY 
2.1. Data Preparation 

In the blind search we analyzed Fermi-LAI data recorded 
between 2008 August 4 and 2014 April 6. The Fermi Sci¬ 
ence Tool^Jwere used to extract Pass 8 source class photons, 
which were analyzed using the P8_SOURCE_V3 instrument 
response functions (IRFs)rJ We used gtselect to select 
photons with reconstructeaairections within an 8° region of 
interest (ROI) around 3FGL J1906.6+0720, photon energies 
> 100 MeV and zenith angles < 100°. We only included pho- 


Hiroshima, Hiroshima 739-8526, Japan 

60 Institucio Catalana de Recerca i Estudis Avangats (ICREA), 
Barcelona, Spain 

61 Department of Physics and Department of Astronomy, University of 
Maryland, College Park, MD 20742, USA 

62 NASA Postdoctoral Program Fellow, USA 

63 Funded by contract FIRB-2012-RBFR12PM1F from the Italian Min¬ 
istry of Education, University and Research (MIUR) 

64 http://tinyurl.com/fermipulsars 

65 http 7/www. einsteinathome . orq 

66 http:77fermi.gsfc.nasa.gov/ssc/data/analysis/ 
software 

D/ The Science Tools, IRFs and diffuse models used here are internal pre¬ 
release versions of the Pass 8 data analysis. Our results did not change sub¬ 
stantially with the final release versions. 













































3 


tons detected when the LAT was working in normal science 
mode, and with rocking angle < 52°. 

To assign photon weights representing the probabil¬ 
ity of each photon having been emitted by the target 
source ( |Kerr| |2011| ), we performed a likelihood spectral 
analysis using the point like package. We built a 
source model by including all 3FGL catalog sources lo¬ 
cated within 13° of 3FGL J1906.6+0720, while allow¬ 
ing the spectral parameters of point sources within 5° to 
vary. We modeled the gamma-ray spectrum of this source 
with an exponentially cu toff pow er lawy typical of gamma- 
ray pulsar spectra ( [Nolan et al.| |2012| . We used the 
template_4years_P8_V2_scaled .fits map cube and 
isotropic_source_4years_P8V3 template to model 
the Galactic diffuse emission and the isotropic diffuse back¬ 
ground respectively]^] The normalization parameters of both 
diffuse components were left free. Finally, the photon weights 
were computed using gtsrcprob, based on the best-fit 
source model resulting from the likelihood analysis. 


2.2. Blind Search Method 

For the blind search, we assumed a canonical isolated pulsar 
model, making it necessary to search in four parameters: spin 
frequency, /, spin-down rate, /, R.A., a and deck, S. 

The basis for most blind searches for gamma-ray pulsars 
is the well-known multistag e scheme based around an initial 
semicoherent search (e.g., Atwood et al . 2006} Pletsch et al. 


2012a). For this survey, we implem ented t he form of the mul 
tistage search scheme described in Pletsch & Clark ( 2014| ), 


where the initial semicoherent stage uses a lag-window of du 
ration 2 21 s « 24 days. 

Notably, this survey incorporates an intermediate semico¬ 
herent refinement step, with a longer (more sensitive) lag- 
window of 2 22 s « 48 days, reducing the parameter space 
around each first-stage candidate to be searched in the final 
fully-coherent follow-up step. This improves the efficiency 
of the follow-up stage, and allows the search to “walk” away 
(in all 4 search parameters) from the original location of the 
candidate if necessary. 

Figure ^illustrates the importance of these new techniques. 
In the blind survey, we searched a conservatively large cir¬ 
cular region around the 3FGL sky location with a radius 50% 
larger than the 3FGL 95% confidence region. As evident from 
Figure]!] the pulsar lies far outside the original source’s confi¬ 
dence region, and also outside our search region. We therefore 
owe its detection to the large resolution of the semicoherent 
step, and the flexibility of the follow-up steps, which allow for 
signals to be detected despite a large offset between the signal 
parameters and the search location. 

The most significant pulsar candidates from the blind 
search were automatic ally refined using the H -test statistic 
(de Jager et al. |1989| . This revealed an interesting candi¬ 
date; however the measured signal-to-noise ratio (S/N) was 
slightly below the detection threshold for a blind search in¬ 
volving a very high number of trials. Upon manual inspection, 
clear pulsations were observed in the photon data after April 
2010; however the phase of these pulsations was not constant, 
and exhibited “wraps” in which the pulsation phase quickly 
jumped by one full rotation. These features indicated that the 
canonical isolated pulsar model used for the blind survey was 


68 http fermi .gsfc.nasa. gov/ssc/data/access/lat%| 
BackgroundModels.html 



Fig. 1.— Sky location of PSR J1906+0722 and positional offset from the 
catalog location. The dotted ellipse shows the 3FGL 95% confidence region. 
The dashed ellipse shows the region in which the search grid (crosses) was 
constructed for the initial semicoherent search stage. The filled area shows 
the region that can be reached by the follow-up stage as it moves away from 
the initial candidate location. The la region from the timing solution, shown 
by the solid ellipse, is highlighted in the inset. 


insufficient, and hence follow-up studies were required to de¬ 
scribe the pulsar’s rotation over the entire dataset. 

3. FOLLOW-UP ANALYSIS 

Before carrying out follow-up analyses, we extended the 
dataset to include photons observed until 2014 October 1 and 
increased the ROI to 10°. To speed up the timing procedure 
computations, we discarded photons with a probability weight 
below 5%. 


3.1. Glitch Identification 

Since pulsations were initially only detected during the fi¬ 
nal 4 years of data, the first step was to identify any glitches 
within the observation time. To achieve this, we searched the 
local {/, /} space around the observed signal, in 150-day seg¬ 
ments with approximately 90% overlap, using the Qm ~test 
( [Bickel et al.|2008t|Pletsch & Clark|2014) , 


Qm = 2M 


EtlhlVn 

E m i 1 2 
n= 1 \ a n\ 


( 1 ) 


where a n and V n are the Fourier coefficients of the measured 
pulse profile and the coherent power at the nth harmonic re¬ 
spectively. Using the Qm test method to weight the contri¬ 
butions from each harmonic, as opposed to the commonly 
used H- test, offers a significant sensitivity improvement ( |Be-| 
ran|1969| ), making it particularly useful when analyzing weak 
pulsar signals. For this step, we included the first 10 Fourier 
coefficients with appreciable power from a segment of the 
data in which the signal was reasonably stable. Using the re¬ 
sults of this scan, shown in Figure [2] an initial ephemeris was 
produced for the timing procedure"described in the following 
section. 
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Fig. 2.— Evolution of the PSR J1906+0722 signal including the glitch at MJD 55067. Left: Phase-time diagram where each point represents one photon, 
with the intensity representing the photon weight. Center and right: The Qio-test (shown by the color-bar) calculated over small ranges in {/,/}, centered on 
the pre-glitch parameters, in overlapping 150-day segments, and maximized over / and / respectively. The dashed line indicates the maximum likelihood timing 
solution. 


3.2. Timing Analysis 

To accurately estimate the pulsar’s rotational, glitch and sky 
locatio n paramete rs we used a variation of the timing method 
used by |Ray etak ( |2011| ), based on unbinned likelihood maxi¬ 
mization. For all N photons in the dataset, with weights {wj}, 
we assigned a rotational phase </> = <f>(tj,u ), determined by 
the photon’s arrival time, tj and the set of model parameters, 
denoted by the vector u. For a template pulse profile, F (</>), 
the likelihood is 


N 

A«) = n \- w i F (^(9> u )) + (! - w i)} • (2) 

i= i 

We first constructed a template pulse profile from the (back¬ 
ground subtracted, see Figure [3]) photons within a sub-section 
of the data set in which the initial ephemeris was believed to 
be accurate. When timing PSR J1906+0722 we used a tem¬ 
plate pulse profile consisting of 3 wrapped Gaussian functions 
( [Abdo et al.[|2013p , which were fit by maximizing the likeli¬ 
hood within the segment. 

With a template profile at hand, we then estimated the pul¬ 
sar’s parameters (given in Table [I]) by varying them around 
their initial estimate to maximize the likelihood over the en¬ 
tire dataset. The result is a likelihood maximization which 
is unbinned in both phase (via the template profile) and time. 
This avoids the need to construct a set of data subsegments for 
pulse times of arrival (TOA) determination. This is especially 


beneficial for faint pulsars, which require longer subsegments 
(and hence fewer TOAs) to ensure the S/N is large enough in 
each for accurate TOA measurement. Subsequently, using the 
most likely parameters, the template profile was updated and 
the process was iterated to maximize the overall likelihood. 

To explore the multi-dimensional parameter space we_used 
the MultiNest nested sampling algorithm (Feroz_et al. 
2013| ), which offers high sampling efficiency, and allows po” 


terior distributions to be calculated as a by-product. 

The timing procedure was carried out in two stages: firstly, 
all timing parameters were allowed to vary. Due to the short¬ 
ness of the pre-glitch segment, the uncertainties in the glitch 
parameters dominated those of the remaining timing param¬ 
eters. We therefore fixed the glitch parameters at their maxi¬ 
mum likelihood values, and fit again for the remaining timing 
parameters. 

When timing radio pulsar glitches, Yu et al. ( |2013| noted 
that unique solutions for glitch epochs could not be found for 
large glitches occurring during an interval between two ra¬ 
dio observations. We observe a similar effect here, although 
our limiting factor is the photon flux. When phase folding, 
a full rotation can be lost/gained if the offset between the 
model glitch epoch and the true glitch epoch is more than 
1/A/ « 0.3 days; however an average of only 1.4 weighted 
photons are observed from the pulsar within this time, mak¬ 
ing this phase wrap simply undetectable. We assumed that no 
phase increment occurred at the glitch, and found that the pos- 
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TABLE 1 

Parameters for PSR J1906+0722 


Parameter 

Value 

Range of Photon Data (MJD) . 

Reference epoch (MJD) . 

54682-56931 

55716 

Timing Parameters 

R.A., a (J2000.0). 

Deck, S (J2000.0) . 

Frequency, / (Hz) . 

1st frequency derivative, /, (Hz s _1 ) ... 

2nd frequency derivative, /, (Hz s -2 ) 

Glitch epoch a (MJD). 

Permanent / glitch increment, A/ (Hz) 

Perm. / glitch increment, A/ (Hz s -1 ) 
Decaying / glitch increment, A/ f ] (Hz) 

Glitch decay time constant, r c \ (days) .. 

19 h 06 m 31 s .20(1) 
+07°22'55"8(4) 
8.9666688432(1) 
-2.884709(2) x 10“ 12 
3.18(1) x 10- 23 

55067+g 

4.033(1) x 10“ 5 
-2.56(3) x 10“ 14 
3.64(9) x 10“ 7 

221 ( 12 ) 

Spectral Properties 

Spectral index, T . 

Cutoff energy, E c (GeV) . 

Photon flux b , F 100 (photons cm -2 s -1 ) 

Energy flux b , G 100 (erg cm -2 s -1 )_ 

1.9 ± 0.1 

5.5 ± 1.2 

(1.1 ± 0.3) x 10“ 7 
(7.3 ± 1.3) x IQ -11 

Derived Properties 

Period, P (ms) . 

1 st period derivative, P(ss _1 ) . 

Weighted 77-test. 

Characteristic age c , r c (kyr) . 

Spin-down power c , E (erg s _1 ). 

Surface P-held strength 0 , Bs (G). 

Light-cylinder P-held 0 , 7 ?lc (G). 

Heuristic distance 0 , d h (kpc). 

111.524136498(1) 
3.587895(2) x HT 14 

731.2 

49.2 

1.02 x 10 36 

2.02 x 10 12 

1.34 x 10 4 

1.91 


Note. — Values for timing parameters are the mean values of the 
marginalized posterior distributions from the timing analysis, with lcr un¬ 
certainties in the final digits quoted in parentheses. 

a Glitch model parameters are defined in Edwards et al. ( 2006), with the 
correction noted by |Yu et ak|(2013f 

b Fluxes above 100 MeV, 7\oo and 6 + 00 , were calculated by extrapolation 
from the E > 200 MeV spectrum. 

c Derived pulsar properties are defined in |Abdo et al.| ( |2013) . The heuristic 
distance, d^ = (Lij;/47rGioo) 1//2 , is calculated from the heuristic luminos¬ 
ity, Lf , described therein. 


terior distribution for the glitch epoch features several bands, 
separated by 1/ A/. Due to the multi-modal shape of the pos¬ 
terior distribution, in Table [Tl we report the glitch epoch that 
results in the maximum likelihood and the 95% credible in¬ 
terval. 

The inclusion of an additional nearby source in the source 
model and raising the energy threshold to 200 MeV when cal¬ 
culating the photon weig hts for PSR J1906+0722 increased 
the S/N (see Section [33] ). Therefore the timing analysis was 
repeated with the updated photon weights, and the results are 
given in Table [T] The time versus rotational phase diagram 
based on this timing solution is shown in Figure [2] and the 
integrated pulse profile is displayed in Figure [3] Through 
these refinement and timing procedures, the initialcandidate’s 
Q 10 -test S/N (Pletsch & Clark] l2014 p was increased from 
0 10 = 6.86 to the highly significant value of #10 = 16.55 
given by the final timing solution. 


_ b9 Th e exp ectation value for Q M under the null hypothesis, Eq [Q m ], in 

|Pletsch & Clark|(|2014) contains an error. S/Ns reported in this work were 
calculated using the corrected value of 4M. 



Fig. 3.— Top panel: weighted pulse profile of PSR J1906+0722 given 
by the timing solution. Lower panels: weighted pulse profiles in increasing 
energy bands. In each p anel, the estimated background level, calculated from 
the photon weights ((Guillemot et al.|2012), is shown by the dashed line . The 
dot-dashed lines mark the oll-pulse phase interval used in Sect ion |3. 3 1 The 
error bars show lcr statistical uncertainties ([Pletsch et al.|2012af Tne pulsed 
fraction, p, and pulsed S/N, # 10 , in each energy band is also shown. 


3.3. Off-pulse Analysis 

Fitting an exponential cutoff model to the spectrum of 
PSR J1906+0722 revealed a relatively high cutoff energy 
compared to typical gamma-ray pulsars ( E c = 6.5 db 0.9 GeV), 
suggesting that the spectrum could be contaminated by the 
presence of a nearby source as was also noted by |Xing &| 
[Wang|pOT4] ). 

To investigate this possibility, we analyzed the off-pulse 
part of the data using photons with energies between 200 MeV 
and 300 GeV. A residual test statistic (TS) map for the off- 
pulse data (see Figure [3]) revealed an excess (0.28 ± 0.02)° 
away from PSR J1906+0722, at (a, 5) = (286.84°, 7.15°), 
with a TS value of 288. 

Modeling this secondary source with a power-law spec¬ 
trum, we added it to the spectral model for the region, keep¬ 
ing its location fixed from the off-pulse analysis, but leaving 
its normalization and spectral index free, and analyzed again 
the full phase interval data. As a result, we found that the 
log-likelihood value increased slightly, and the new photon 
weights increased the S/N of the pulsations from #10 = 16.38 
to #10 = 16.55. 

The low energy threshold of 200 MeV was chosen to pro¬ 
vide improved angular resolution in order to better separate 
the pulsar emission from that of the new source. When lower 
energy (100-200 MeV) photons were included in the spec¬ 
tral analysis, the pulsation S/N calculated with the resulting 
photon weights decreased, suggesting that source confusion 

































































6 


Clark et al. 


at low energies leads to a less reliable source model. 

Figure [4] shows TS maps and spectral energy distributions 
for PSR J1906+0722 and the new source found in this off- 
pulse analysis. The integrated energy flux of the secondary 
source above 100 MeV is 4.341 q!67 x 10 -11 erg cm -2 s -1 
with a spectral index of 2.17 d= 0.07. 

The best-fitting location of the secondary source is very 
close to the western edge of the supernova remnant (SNR), 
G41.1-0.3 (3C 397,|Safi-Harb et al.|2005). |Jia ng et al.|p0l0| 
observed a molecular cloud interacting with the SNR at this 
location; it is possible that we are observing gamma-ray emis¬ 
sion resulting from this interaction. 


4. ANALYSIS IN OTHER WAVELENGTHS 
4.1. Radio and X-ray Observations 

In probing for radio emission from PSR J1906+0722, we 
carried out a 120-minute follow-up observation with the L- 
band (1.4 GHz) single-pixel receiver mounted on the 100 m 
Effelsberg Radio Telescope in Germany. The gamma-ray- 
timing ephemeris allowed us to search the data over dis¬ 
persion measure (DM) only. No evidence for radio pulsa¬ 
tions was found. Assuming a 10% pulse width, bandwidth 
A F = 150 MHz, telescope gain G = 1.55, n p = 2 po¬ 
larization channels, system temperature T sys = 24 K, digi¬ 
tization factor /3 = 1.2 and a signal-to-noise threshold of 
5, by the radio meter equation (Equation (A 1.22), Lorimer &| 


Kramer 2005 


p.265), we computed a flux density limit of 
21 /iJy. While this is below_the conventional radio-quiet 
level of 30 //Jy ( |Abdo et~aL|[2013| ), we note that the nearby 
LAT-discovered pulsar PSR J1907+0602 has been ob served 
in radio observations with a flux density of just 3.4 /iJy ( |Abdo| 
|et al.|2010l ), and would therefore not have been detected in 
this radio search. 

To check for a possible X-ray counterpart, we an alyzed 
a 10 ks observati on with Swiff s X-ray Telescope (Stroh 


& Falcone 2013| ). No counterpart source was detected, 


with an unabsorbed-flux (0.5-10 keV) upper limit of 2 x 


10 15 erg cm 1 s 1 * * at the pulsar position. This limit yields 
a gamma-ray-to-X-ray flux ratio of > 365, or an efficiency 
Lx/E < 8. 7xl0~ 5 at distance dh, similar to othe r gamma- 
ray pulsars ( [Marelli et al.|2011[|Pletsch et al.|2012a| ). 


4.2. Possible SNR Associations 

There are 4 known SNRs lying within l^_from the tim¬ 
ing position of PSR J1906+0722 (|Green||2014|. There is 
strong evidence that the closest of these, G41.1 —0 .3, is a Type 
la SNR fro m a Chandrasekhar mass progenitor (Yamaguchi 
|et al. |2015] ), making it unlikely to be the birthplace of a pulsar. 
Each of the remaining nearby SNRs lies closer to other young 
pulsars than to PSR J1906+0722 (G41.5+0.4 and G42.0-0.1 
to PSR J1906+0746; G40.5-0.5 to PSR J1907+0602), mak¬ 
ing a physical association between any of these difficult to 
verify. Kick-velocity requirements based on the pulsar’s char¬ 
acteristic age and heuristic distance do not rule out any of 
these SNRs as the birthplace of the pulsar. 


5. DISCUSSION 

Despite several years of attempts, the identification of 
2FGL J1906.5+0720 remained elusive. Now that this source 
has been identified as PSR J1906+0722, we here investigate 
potential reasons for the failure of previous searches to detect 
it. 


Perhaps the most significant source of difficulty in the de¬ 
tection of PSR J1906+0722 was the large positional offset 
between its 3FGL catalog position and its true position. This 
offset, which could only be acc omm odated by the new follow¬ 
up method outlined in Section |2.2| is most likely due to the 


3.3 


presence of the secondary source described in Section 

The close proximity of PSR J1906+0722 to the Galactic 
plane (b = 0.03°) likely also hindered its detection, as the 
large majority of the weighted photons can be attributed to 
the background. From the pulse profile shown in Figure [3j 
we estimate that the pul sed fraction of the total weighted pho¬ 
ton flux (as defined in Plet sch & Clark| ( |2014| )) is as low as 
6%. This low pulsed fraction leads to a low observable S/N, 
making detection more challenging. 

A further complication for detecting PSR J1906+0722 was 
the presence of the glitch about one year into the Fermi 
mission. This glitch is among the largest detected from a 
gamma-ray pulsar in terms of relative magnitude (A/// « 
4.5 x 10 -6 ) (Pletsch et al. 2012b). In previous searches us¬ 
ing a shorter total observation time, the data segment after the 
glitch represented a much shorter fraction of the total obser¬ 
vation time. As the time interval covered by Fermi's observa¬ 
tions since 2008 August 1 continues to increase, the existence 
of a long timespan in which a pulsar’s signal is stable becomes 
ever more likely. The increase in the weighted photon flux 
offered by the Pass 8 analysis ( [Atwood et al.12013] ) further in¬ 
creases the observable S/N throughout the observation time, 
and results in searches that are not only more sensitive overall 
( [Laffon et al.|2015| , but also more robust against glitching or 
noisy pulsars. 

The ability to detect young gamma-ray pulsars in blind 
searches can be of significant importance to the overall study 
of energetic pulsars. For example, Ravi et al. ( [2010 ) use the 
observed population of radio-quiet pulsars to investigate the 
dependence of properties of pulsar emission geometries on 
the spin-down energy, E. Since pulsars with a high E tend to 
exhibit timing noise and glitches (which do not typically af¬ 
fect radio searches), they are hard to find in gamma-ray data, 
where long integration times are required. Advanced search 
methods that can detect complicated signals such as that from 
PSR J1906+0722 are therefore crucial for reducing a poten¬ 
tial bias against young, energetic and glitching pulsars in the 
radio-quiet pop ulation. As noted by Ab do et al.| ( [2013 ) and 


Caraveo ( 2014), such pulsars are indeed lacking in the Fermi 
pulsar sample. 
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Fig. 4.— Top panels: test statistic (TS) sky-maps of the PSR J1906+0722 region above 200 MeV. Each pixel shows the TS value for a point source located at 
the pixel position. The cross represents the timing position of PSR J1906+0722. The light ellipse shows the 95% confidence region of the 3FGL source, the blue 
ellipse shows the 95% confidence region of the new secondary source, and the darker ellipses show the approximate extents o f ne arby SNRs. Bottom panels: 
Spectral energy distributions for the full-pulse interval. The solid curves present the results of the likelihood analyses of Section |T3| 
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